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The existence of a second quadratic inviscid invariant, the helicity, in a turbulent flow leads to 
coexisting cascades of energy and helicity. An equivalent of the four-fifth law for the longitudinal 
third order structure function, which is derived from energy conservation, is easily derived from 
helicity conservation pIM . The ratio of dissipation of helicity to dissipation of energy is proportional 
HN ' to the wave-number leading to a different Kolmogorov scale for helicity than for energy. The 

^^ , Kolmogorov scale for helicity is always larger than the Kolmogorov scale for energy so in the high 

^^ ■ Reynolds number limit the flow will always be helicity free in the small scales, much in the same 

way as the flow will be isotropic and homogeneous in the small scales. A consequence is that a pure 
helicity cascade is not possible. The idea is illustrated in a shell model of turbulence. 
Few exact results regarding fully developed turbulence have yet been derived. The most celebrated being Kol- 
mogorovs four-fifth law j3|]. The four-fifth law is based on the fact that energy, which is an inviscid invariant 
of the flow, is transferred through the inertial range from the integral scale to the dissipation scale. The four- 
fifth law, (8v(l)?,) = —(4/5)el, states that the third order correlator associated with energy flux equals the mean 
energy dissipation. As noted recently HB] in the case of helical flow a similar relation exists for the transfer 
J^T of helicity leading to an other scaling relation for a third order correlator associated with the flux of helicity, 
t—( ' (<5v||(Z) • [v_l(j") x v±(r + I)]) = (2/15)Sl 2 , where 8 is the mean dissipation of helicity. This relation is called the 
f^ . 'two-fifteenth law' due to the numerical prefactor. This establishes another non-trivial scaling relation for velocity 
< ^> ' differences in a turbulent helical flow. 



The coexistence of cascades of energy and enstrophy is prohibited for high Reynolds number flow in 2D turbulence. 
The reason for this is that the enstrophy dominates at small scales such that the ratio of energy - to enstrophy 
dissipation vanishes for high Reynolds number flow. The Kolmogorov scale k^ for enstrophy dissipation is determined 

from the energy spectrum E(k) ~ fc -3 and the kinematic viscosity v by C, — v J z dkk 4 E(k) ~ vk 2 z => kz ~ v~ x l 2 . 



The energy dissipation is e = v J z dkk 2 E(k) ~ flog kz ~ — (l/2)i/logi/ — > for u — ► 0. Consequently energy is 
cascaded upscale in 2D turbulence. 

The situation in 3D turbulence is different. Here coexisting cascades of energy and helicity are possible p|. However, 
-— . the same type of dimensional argument as for the cascades of energy and enstrophy in 2D turbulence applies. The 
helicity density is h — u%Wi, where Wj = £ijkdjUk is the vorticity. The mean dissipation of helicity is Dh — v{djUidjUi). 
Disregarding signs this can spectrally be represented as 



H ! D H ~v j dkk 3 E(k) ~ vk 7 J 3 ~ f" 3/4 , (1) 

where fc^ 1 = r\ is the Kolmogorov scale and we have used E(k) ~ fc~ 5 / 3 and Ue ~ f~ 3 / 4 . This means that for 
high Reynolds numbers flow the dissipation of helicity will grow as i?e 3 / 4 . Since the mean dissipations of energy e 
and helicity 8 are determined by the integral scale forcing the growth of helicity dissipation with Reynolds number 
is apparently in conflict with the assumption of a constant energy dissipation in the limit of vanishing viscosity. 
This is not a true problem because helicity is non-positive, and the viscous term in the equation for the helicity 
v(uidjjU>i + LOidjjUi) can have either sign. So in the high Reynolds number limit there must either be a detailed 
balance between dissipation of positive and negative helicity or the energy cascade is blocked jq]. In the rather 
artificial case of a shell model where only one sign of helicity is dissipated by hyper-viscosity, the energy cascade is 
indeed prevented all together similar to the case of forward energy cascade of energy in 2D turbulence jq] . 

In a helical flow (S ^ 0) the dissipation of helicity defines a scale different from the Kolmogorov scale r]. This we 
will call the Kolmogorov scale £ for dissipation of helicity. 

Following K41, the Kolmogorov scale r\ for energy dissipation is obtained from e ~ Su^/7] ~ v8u 2 /ri 2 => 77 ~ 
(f 3 /^) 1 / 4 , where 6ui is a typical variation of the velocity over a scale I. The Kolmogorov scale £ for dissipation 
of helicity is defined as the scale where the helicity dissipation is of same order as the spectral helicity flux. With 
dimensional counting we have 8 ~ v8u\/£, 2 and using 8ui ~ (Is) 1 ' 3 we obtain 



e-o^ 2 /* 3 ) 177 - (2) 

Now it is clear why (P leads to a wrong conclusion for the mean dissipation of the helicity S. The integral will not 
be dominated by contributions from fcg but contributions from ku = l/£, 

D H =~S~ vk 1 ^ ^k H ~ v- 3 ' 7 . (3) 

The ratio of the two Kolmogorov scales is then (r)/(,) = (fc^/fes) ~ jy~ 3 / 7 + 3 / 4 — ^ 9 / 28 _» o for v — > 0. Thus for high 
Reynolds number helical flow the small scales will always be non-helical and a pure helicity cascade is not possible. 

On the other hand for scales I < £ the ratio of dissipation of energy and helicity is proportional to I, De/Dh ~ I, 
which means that helicity dissipation dominates and the dissipation of positive and negative helicity must balance. 
The reason for the flow to be non-helical on small scales is different from the reason why the flow tends to be isotropic 
on small scales even though the integral scale is non-isotropic. The reason for the small scales to be isotropic is that 
the structure functions associated with the non-isotropic sectors scale with scaling exponents that are larger than 
those of the isotropic sector and thus becomes sub-leading for the flow at small scales independent of the dissipation 

§• 

The physical picture for fully developed helical turbulence is then that S and e are solely determined by the forcing 
in the integral scale. There will then be an inertial range with coexisting cascades of energy and helicity with third 
order structure functions determined by the four-fifth - and the two-fifteenth laws. This is followed by an inertial 
range between £ and r\ corresponding to non-helical turbulence, where the dissipation of positive and negative helicity 
vortices balance and the two-fifteenth law is not applicable. 

In order to test these ideas in a model system we investigate the role of helicity and the structure of the helicity 
transfer in a shell model. 

Shell models are toy-models of turbulence which by construction have second order inviscid invariants similar to 
energy and helicity in 3D turbulence. Shell models can be investigated numerically for high Reynolds numbers, in 
contrast to the Navier-Stokes equation, and high order statistics and anomalous scaling exponents are easily accessible. 
Shell models lack any spatial structures so we stress that only certain aspects of the turbulent cascades have meaningful 
analogies in the shell models. This should especially be kept in mind when studying helicity which is intimately linked 
to spatial structures, and the dissipation of helicity to reconnection of vortex tubes H . So the following only concerns 
the spectral aspects of the helicity and energy cascades. 

The most well studied shell model, the GOY model M, is defined from the governing equation, 

w'n = ik n {u n+2 u n+ i - -u„+iu„-i + —— Un __ lUn _ 2 )* - uk 2 n u n + f n (4) 

with n = 1, ...,N where the u rl 's are the complex shell velocities. The wavenumbers are defined as k n = A™, where A 
is the shell spacing. The second and third terms are dissipation and forcing. The model has two inviscid invariants, 
energy, E — J2 n ^n — J2 n \u n \ 2 , and 'helicity', H = J2 n Hn = Sn( e ~ l)~"l u n| 2 - The model has two free parameters, 
A and e. The 'helicity' only has the correct dimension of helicity if |e — l|~ n = k n => 1/(1 — e) — A. In this work we 
use the standard parameters (e, A) = (1/2, 2) for the GOY model. 

A natural way to define the structure functions of moment p is through the transfer rates of the inviscid invariants, 

s*(k n ) = ((nf y'^k-^ (5) 

s?(kn) = ((n^/ 3 >fc- 2 */ 3 (6) 

The energy flux is defined in the usual way as II„ = d/dt| n .;.(X)m=i ^m) where d/dt\ n .i. is the time rate of change 
due to the non-linear term in (0) . The helicity flux 11^ is defined similarly. By a simple algebra we have the following 
expression for the fluxes, 

<n r f) = (l-e)A„ + A„ +1 =e (7) 

(I#) = (-l) n k n (A n+l - A„) = S (8) 

where A„ = fc n _iim(u„_iw n u„ + i), e and 5 are the mean dissipations of energy and helicity respectively. The first 
equalities hold without averaging as well. These equations are the shell model equivalents of the four-fifth - and the 
two-fifteenth law. Kadanoff et al. defined a third order structure function as, 

S 3 = Im{u n -xu n u n +\) = A„/fc„-i (9) 



to avoid the spurious (specific to the GOY model) period 3 oscillation. Using this we obtain from ([7]) and (g) a scaling 



relation for Si 
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" (1 - e/2)k n 



(e-(-l) n 5/k n ). 



(10) 



The last term in the parenthesis is sub-leading with period two oscillations. When S — the sub-leading term 
disappears and the scaling from the four-fifth law is obtained, figure 1. The relation (0) is the scaling relation 
corresponding to the sub- leading term, which survives due to detailed cancellations between the two terms A„+i and 
A n of the leading term corresponding to (m). The case e = and 5 ^ would, aside from the period two oscillation, 



k 



-2/3 



However, this 



correspond to a helicity cascade with the scaling obtained from dimensional counting u T , 
situation is, as we will show shortly, not realizable. 

The mean dissipations e and S are from energy and helicity conservations identical to the mean energy and helicity 
inputs which from M) are, 



e = ^(/„<)+c.c. 



(11) 



and 



s^y^(-i) n k n (f n ui)+c.c. 



(12) 



so e and 8 are not independent. The forcing can be chosen in many ways. A natural choice is /„ = /„/«*, where 
/° is independent on the shell velocities. Then we have, e = J2 n <n fn an d ^ = J2 n <n ( — 1)"^™/°! n i indicates the 
end of the integral scale. By choosing the coefficients, stochastic or deterministic functions of time, this last sum 
can vanish identically, which is referred to as helicity free forcing. The simulations shown in figure 1 are performed 
with the forcing, /£ = 1CT 2 (1 + i) and f% = -Af£/\ with A = and A = l, corresponding to (e,6) = (0.01, 0) and 
(e ,6) = (0.01,0.08) respectively. 

Helicity is not positive and is dissipated with opposite signs for odd and even shells. If we consider the third order 
structure function associated with the helicity transfer as defined by (H) we see (figure 2) period two oscillations 
growing with n. This period two oscillation is due to the dissipation and not the non-linear transfer. The helicity flux 



(nZ) = 8-(D n ), 

where D n is the helicity dissipation at shells m < n: 

n 
Dn = Y, K-l)™fc 



3 l U I 2 

U/rn. 



(13) 



(14) 



In the inertial range for energy transfer we have the Kolmogorov scaling u r , 
estimated, 



D„ 



V (-irkW~\v* 



(_l)n A 7n/3 . 
A 7 /3 + 1 



1 



.-1/3 



so the helicity dissipation can be 



(-1)»# 8 . 



(15) 



Figure 3 shows |(n^)| and (11^) as functions of wave number. The scaling (fL5|) of the helicity dissipation is the straight 
line, the horizontal dashed line is 6. The inertial range for helicity transfer is to the left of the crossing of the two lines. 
The crossing is the Kolmogorov scale for helicity transfer Kh, which does not coincide with the Kolmogorov scale 
for energy transfer, Ke- The 'pile- up' for k larger than Kh was earlier interpreted as a bottleneck effect pQ]. It is a 
balance between positive and negative helicity dissipation. The forcing f n = /„/u„ can potentially cause numerical 
trouble when \u n \ becomes small. It is easy to see that the linear equation for (real) shell velocity u n , neglecting 
the non-linear transfer, ii n = f /u n will create a finite time singularity. This is not the case for the forcing suggested 
by Olla at two shells, /„ = (aE n+1 u n ) / (E n + E n+1 ) and f n+1 = (bE n u n +i) / (E n + E n+X ), where a and b are 
constants determining the ratio of energy - to helicity input. The coupled set of equations, (u n = f n , ii n +i = fn+i) 
is integrable (solve for y = u n /u n+ i), and has no finite time singularities. Using this forcing we performed a set of 



simulations with constant energy input e = 0.01 and varying helicity input 5 = (0.0001,0.001,0.005,0.01,0.08). In 
figure 4 the spectra of the absolute value of the helicity transfer normalized with S are plotted against wave number 
normalized with Kh ■ Kh is in each case calculated from (||), and a clear data collapse is seen. 

In summary the simulations with the GOY shell model suggest a new Kolmogorov scale for helicity, always smaller 
than the Kolmogorov scale for energy. Thus there exist two inertial ranges in helical turbulence, a range smaller than 
Kh with coexisting cascades of energy and helicity where both the four-fifth - and the two- fifteenth law applies, and 
a range between Kh and Ke where the flow is non-helical and only the four-fifth law applies. 
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FIGURE CAPTIONS 

Fig. 1 The third order structure function S^ as calculated from (|jj) in the cases 6 > (crosses) and 5 — (diamonds). 
In the case of helicity free forcing the modulus 2 oscillations disappears. In the two runs we have 25 shells, 
v = 1(T 9 , /„ = 0.01(1 + i)(6 n , 2 /u2 - A6 n ,a/2u$) with A = 0, 1 respectively. 

Fig. 2 The helicity flux (11^) in the case 5 > 0. The same curve is multiplied by 1000 and over-plotted in order to see 
the inertial range. The period 2 oscillations in the helicity transfer comes from the helicity dissipation. 

Fig. 3 The absolute values of the helicity flux |(II^)| (diamonds) show a crossover from the inertial range for helicity 
to the range where the helicity is dissipated. The line has a slope of 7/3 indicating the helicity dissipation. The 
dashed lines indicate the helicity input <5. The crosses is the helicity flux in the case 3 = where there is no 
inertial range and Kh coincides with the integral scale. The triangles are the energy flux (IF^). 

Fig. 4 Five simulations with constant viscosity v = 10~ 9 , constant energy input e = 0.01 and varying helicity input 
1> = (0.0001,0.001,0.005,0.01,0.08) are shown. The absolut values of the helicity flux \(U^)\ divided by 5 is 

plotted against the wave number divided by Kh = (^ 3 e 2 /5 )~ 1 ^ 7 , which is obtained from (0) neglecting O(l) 
constants. A clear data collapse is seen. 
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